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Abstract 

We consider the Landau-Ginzburg- Wilson Hamiltonian with 0(n) x 0(m) symme- 
try and compute the critical exponents at all fixed points to 0{n~ 2 ) and to 0(e 3 ) in a 
e = 4 — d expansion. We also consider the corresponding non-linear cr-model and de- 
termine the fixed points and the critical exponents to 0(e 2 ) in the e = d — 2 expansion. 
Using these results, we draw quite general conclusions on the fixed-point structure of 
models with 0(n) x 0(m) symmetry for n large and all 2 < d < 4. 
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1 Introduction. 



The critical behavior of frustrated XY and Heisenberg spin systems with noncollinear or- 
der has been the subject of many recent theoretical studies, where the standard tools of 
renormalization-group (RG) theory have been applied to field theories which were conjec- 
tured to be appropriate for the description of the systems under investigation (see, e.g., 
Refs. fl],|2| for reviews on this issue). 

The critical behavior of these systems is rather controversial. Indeed, while experimen- 
tally there is good evidence of a second-order phase transition^ belonging to a new (chiral) 
universality class, theoretically the issue is still debated. On one side, field-theoretical studies 
based in approximate solutions of the RG equations (ERG) do not find any stable fixed point 
and favor a first-order phase transition [0-0. On the other hand, perturbative field theory 
gives the opposite answer: A stable fixed point is identified with exponents in agreement 



with the experiments [ 10| . Monte Carlo simulations |TT]-[y| do not help clarifying the issue. 
While simulations of the antiferromagnetic 0(n) model on a stacked triangular lattice find 
a second-order phase transition with exponents reasonably near to the experimental ones, 
modified spin systems which supposedly belong to the same universality class apparently 
favor a first-order transition ||13|| f| Note that the existence of a new chiral universality class 
does not exclude the possibility that some systems undergo a first-order transition. Indeed, 
they may lie outside the attraction domain of the stable fixed point and thus belong to 
runaway RG trajectories. In this first-order transition is expected. 

The field-theoretical studies have been focusing either on the so-called Landau-Ginzburg- 
Wilson (LGW) Hamiltonian with 0(n) x 0(2) symmetry or on the corresponding nonlinear 
sigma (NLa) model. In this paper we will study a generalization of these theories, by 
considering general 0(n) x 0(m) Hamiltonians and we will try to understand the nature of 
the fixed points of the theory. In particular, we will relate the LGW and the NLa descriptions 
showing explicitly that the stable fixed points of the two models are exactly the same, as 
conjectured in Refs. [P^TB| , for values of m and n consistent with the existence of a second- 
order phase transition. Moreover, we will clarify the nature of the unstable fixed points. 
However, this analysis will only be valid in the large-n region, where, by using the large- 
n expansion, we will be able to identify nonperturbatively all fixed points of the different 
Hamiltonians. 

Since the large-n expansion plays a major role in our discussion, our results will only be 
valid for n > n(m, d), i.e. in the region of large-n analyticity. At least near four dimensions, 
such a function is conjectured to be identified with the line n + (m, d) on which the LGW chiral 

1 The experimental results are reviewed in, e.g., Refs. 0,0). Essentially, experiments with hexagonal 
peroskvites find a clear second-order phase transition except for CsCuCl3. The results for helimagnetic rare 
earths are instead less clear. We also mention Ref. H where it is shown experimentally that chiral order and 
spin order occur simultaneously, thereby supporting Kawamura's conjecture that chiral transitions are 
different from the standard 0(n) transitions. 

2 Note that not all modified models show a first-order phase transitions. Some of them have a behavior 
that is compatible with a second-order phase transition. However, the measured exponents do not satisfy 
the condition r\ > 0, which must be satisfied in unitary (reflection-positive) models as these are, so that, 
the measured exponents can only be effective ones. This is interpreted as a signal of a first-order phase 
transition. However, there is also the possibility that the results are strongly biased by corrections to scaling 
induced by the constraint. 
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and antichiral fixed points merge. The function n + (m,d) has been the object of extensive 
studies that tried to understand whether, in the physical case d = 3 and m = 2, n(m, d) was 
smaller or larger than three. In this case, studies using various approaches gave n(2, 3) ~ 4 
(Refs. P, |T7| , |I%|] ), 5 (Ref. |J), ps 6 (Ref. [PI). Here, we will provide another determination, 



together with generalizations for other values of m, that substantially confirms previous 
findings, i.e. n(2, 3) ~ 5. Since the results that we will present are essentially adiabatic 
moving from the large-n and small e = 4 — d region, they are not expected to provide the 
(essentially nonperturbative) features of the models in the region below n(m,d). Therefore, 
the fact that n(2, 3) > 3 does not necessarily imply an inconsistency with the field-theoretical 
results of Ref. [fLOfl , where a rather robust evidence for stable chiral fixed points was found 
for 0(n) x 0(2) models with n = 2, 3 in fixed dimension d = 3. Such fixed points are not 
analytically connected with the large-n and small-e criticalities discussed above. 

We will also show that for m > 2 the identification n(m,d) = n + (m,d) may not be 
correct for d near two dimensions. Indeed, in this case a new critical line appears which 
corresponds to the merging of the chiral fixed point with the NL<r antichiral fixed point. 

In order to obtain quantitative predictions for all m and n, we have extended the e = 4 — d 
expansion of the LGW theory to order e 3 and the e = d — 2 expansion of the NLcr model to 
order e 2 for all n and m. Also, we present 0(l/n 2 ) results for the LGW theory. 

The paper is organized as follows. In Sec. ^| we define the general class of models with 
0(n) x 0{m) symmetry that will be considered in the paper and find a general represen- 
tation that is the starting point of the large-n expansion. In Sec. |3] we compute the 0(e 3 ) 
contributions to the critical exponents r\ and v~ x within the e = 4 — d expansion of the LGW 
Hamiltonian. In Sec. |4] we analyze in detail the 1/n expansion of the LGW Hamiltonian 
with 0(ri) x 0(m) symmetry to 0(l/n 2 ), thereby extending the results of Ref. |TJj]. Inter- 
estingly enough, we can determine the large-n expansion of the exponents at all fixed points 
and show explicitly their different physical nature: at the stable fixed point both tensor and 
scalar excitations propagate, while at each unstable fixed point one of the degrees of freedom 
is suppressed. At the Heisenberg fixed point there are only scalar excitations, while at the 
antichiral one, there are only tensor excitations. In Sec. [5] we discuss the 1/n expansion of 
a more general theory in which the coupling to a (gauge) vector field is included, extending 
the results of Ref. [19[]. In Sec. |6] we extend to arbitrary values of m and to C^e 2 ) the e- 
expansion of NLa models, evaluating the unstable fixed point and the coalescence value of n 
under which the two fixed points actually disappear. We also identify the "gauge" criticality 
of the models. In Sec. |7| we draw some general conclusions and present a new determination 
of the function n(m, d). 



2 Models. 

We will consider a non-Abelian gauge model coupled to a scalar field with gauge symmetry 
0(m) and global symmetry 0(n). In particular, we consider a set of m n-dimensional vectors 
4> a = {'Paa}, & = 1, • • • , m , a = 1, ■ ■ ■ ,n, a vector field A"f antisymmetric in a and /3, and 
the Hamiltonian density 

n = \ E(d^ a + 9oAf^) 2 + \r, fe 
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E [Ma ■ - <I>14>1] + \f% + I E A? A?, (2-1) 

a/3 a/3 

where F^ v is the non-Abelian field strength associated with the fields A"". This Hamiltonian 
is gauge-invariant (with local 0(m) invariance) for t Q = 0, and in this case it has already 
been studied^ in Ref. |19|| . For to — and A°f = we obtain a generic LGW Hamiltonian 
density with global 0{m) x 0{n) invariance: 

n = \ E(W 2 + \r Q X> 2 + ±«o (e ^) 2 

+^oE[(^-^) 2 -^]- (2-2) 

4 - a/3 

Assuming n > m, stability requires uq > and wo = wo + (1 — M)vq/N > 0, where 
A/" = min (m, n). 

Other particular cases of the Hamiltonian Q2.1|) are interesting. If we set Wo = t>o and 
r = —v rji/6 and take the limit v — > +oo keeping 771 fixed, we obtain an 0(n) x O(m) 
cr-model coupled to an 0(m) vector field. The Hamiltonian density is given by 

U = \ E(^« + 9oAfd>,) 2 + -A + I E AfAf, (2.3) 

where the fields <f> satisfy the constraint 

4> a -<t>p = Vl8a/3- (2.4) 



This limit is well defined only if n > m, otherwise the constraint ( |2.4[ ) cannot be satisfied. 
In the absence of the kinetic term for the vector field, the Hamiltonian density depends 
quadratically on the vector field that can then be eliminated by integration. We obtain a 
new Hamiltonian of the form 

n = ~ E( - R(n 2?*, t s Y,(<t>«We - WM 2 - (2-5) 



This is the cr-model studied in Refs. [|T^-[TB|. In order to recover the notationsQof Ref. [W 
we set (f) a = ^/rj{e a and r/ 2 = 2t r] 1 /(g^rj 1 + 1 ). Then 

1 m fl \ m 

ff=^iE <^ e « ' d v e « + ( nV2 - Vi ) E ( e « ' d n e pfi (2-6) 

1 a=l ^ Z J a>0 



3 Hikami's couplings labelled by the subscript H, are related to ours by the correspondence: pn = 
(uo - «o)/3, v H = v /6, X H = Ph + 2v H = ito/3. 

4 Notice that our rji are consistent with the couplings employed in Sec. 4 of Ref. but they are twice as 
big as the couplings defined in Appendix B of the same paper, due to a slight inconsistency in the notation 
adopted by these authors. The couplings r\i are related to those of Ref. |14| by r\\ = 1/T + 1/T", 772 = 2/T. 
For easy comparison, the reader should keep in mind that, due to the constraints, e a ■ d^ep = — ep ■ d^e a . 
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where the fields e a are m n-component vectors (or equivalently mxn matrices) with n > m 
subject to the nonlinear constraint 

e a ■ ep = 8 a p. (2.7) 

In order to study the large-n behavior of these models we rewrite the general Hamiltonian 
(p.l|) in such a way that the dependence on the field <p a is quadratic. This is obtained by 
introducing two auxiliary fields: a scalar field S and a symmetric and traceless tensor field 
T a/3 , i.e. such that T Q/3 = T@ a , T aa = 0. By means of these auxiliary fields we can rewrite 
the Hamiltonian ( |2.1| ) as 

H = H eS - — T 2 - ^S 2 + —A 2 + -F 2 } (2.8) 

2 2 2 4 y J 

where 

H eS = W<f> a - xaP ^ ( 2 -9) 

and 

X a ? = -d,d,5 aP + r 5 a(S - 2g Afd, + v T a P + w S5 a/3 + g 2 A^A a ^. (2.10) 

Note that the effective action for the <fi fields is the most general one which is 0(m) covariant. 
Therefore, the analysis of this class of models provides the critical behavior of the most 
general 0(m) x 0(n) theory. 



3 e-expansion for the Landau- Ginzburg- Wilson model 



In this Section we study the LGW Hamiltonian (|2.2|) and report our results for the critical 
exponents and the /^-function to order e 3 , thereby extending the results of Ref. |20| to three 
loops. We consider the massless theory and renormalize it using the MS scheme. We set 



= [Z^u,v)} l ' 2 cf> R , 
u = p e Z u (u,v)N d 1 , 
v = (SZvfav)^ 1 , 



(3.1) 
(3.2) 
(3.3) 



where the renormalization constants are normalized so that Z^(u,v) « 1, Z u (u, v) ~ u, and 
Z v (u, v) « v at tree level. Here Nd is a <i-dependent constant given by iV^ 1 = 2 d - 1 n d / 2 T(d/2). 
Moreover, we introduce a mass renormalization constant Z t (u,v) by requiring Z t T^ 1,2 ^ to be 
finite when expressed in terms of u and v. Here I^ 1 ' 2 ) is the two-point function with an 
insertion of <fi 2 . Once the renormalization constants are determined, we compute the f3 
functions from 



fi u (u,v) = n 
and the critical exponents from 



du 
d/j, 



UQ,V 



dv 
d/j, 



u ,v 



d log Z^ 
dlogfi 



(3.4) 



(3.5) 



UQ,Vo 



4 



Vt 



d\ogZ t 
d log /i 

(2 - V ~ Vt) 



u ,v 

-1 



For the /3-functions we obtains 



-eu 



mn + 8 
6~ 



-u — 



(m — l)(n — 1) / v 

v [u 

3 V 2 



3mn + 14 
12 



-u 



/ll 2 13 5 

+ m - 1 n - 1 d — u uv -\ 1 

v A ; U8 24 36 

33m 2 n 2 + 922mn + 2960 + C(3)(480mn + 2112) 

(m - l)(n - 1) {-4[79mn + 1318 + 768C(3)]w 3 

+ [555wm - 460(m + n) + 6836 + 4032C(3)]m 2 i; 
-2[213mn - 358(m + n) + 1933 + 960C(3)]m; 2 
+[121mn - 309(m + n) + 817 + 216C(3)]w 3 } , 



u 

1728 

v 

3456' 



(3.6) 
(3.7) 



(3i 



m + n — 8 2 5mn + 82 2 
— u v 



-ev + 2uv , 

6 36 
5mn — ll(m + n) + 53 2 13mn — 35(m + n) + 99 3 



-uv 



72 



3456 



|52m 2 ra 2 — 57mn(m + n) — 220Qmn — lll(m 2 + n 2 ) + 4291 (m + n) 



+ [-1416mn + 3216(m + n) - 7392]C(3)} 

\ —39m 2 n 2 + 35mn(m + n) + 1302mn + 36(m 2 + n 2 ) 
864 I v ' v 

+ [768mn - 1824(m + n) + 4896]C(3)} 
yj2g ( ' § m n — 35mn(m + n) — 2114mn + 3182(m + n) 
+ [-1152mn + 2304(m + n) - 10368]C(3)} 



v 3 u 



U 3 V 

"864 



-13mV + 3Q8mn + 3284 + (192mn + 2688)C(3) 



For the critical exponents we obtain: 
mn + 2 



V 



72 



u 2 + (m — l)(n — l)v ( — — | 

V A ; V48 36/ 



(mn + 2)(mn + 
1728 



8084 



2401(m + n) + 5725 



12520 



(3.9) 



We mention that for the particular case n = m = 2 one may derive the corresponding four-loop e- 
expansion from the series reported in Ref. pl| for the so-called tetragonal model. Indeed an exact mapping 
JfTj exists bringing from the LGW Hamiltonian (2.2) with m — n — 2 to the tetragonal model considered in 
Ref. |Q . Note that our renormalized couplings differ from those defined for m — 2 in Ref. || . Kawamura's 
couplings uk, vk are related to ours by u = X2N^uk, v = QNdVK- 
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( m — 1 ) ( — 1 ) r 

v v 2 {2mn — 5m — 5n + 26) 

3456 L 

—6uv(mn — m — n + 10) + 6u 2 (mn + 8) 



(3.10) 



1 mn + 2 (m — l)(n — 1) 5(mn + 2) 9 

- = 2 u + - - -v + — V 

v 6 6 72 



+5(m - l)(n - lW 

v A ; V48 36/ 

(m — l)(n — 1 



u (mn + 2)(hmn + 37) 3 



288 

3v 2 (7mn — 16m — 16n + 79) 



1152 

-w(61mn - 58m - 58n + 550) + 12u 2 (5mn + 37)1 . (3.11) 

As discussed at length in Refs. Hlj|6|,^0l, the critical behavior of these systems depends on 
the values of n and m. In general, the /3-functions admit four solutions: the Gaussian fixed 
point (u* = v* = 0), the 0(mn) Heisenberg fixed point (v* = 0) and two new fixed points 
with nontrivial values of u* and v*, the chiral and the antichiral fixed points. These two 
additional fixed points do not exist for all n and m, but, at fixed m, only for n > n + (m) and 
n < n~(m). The functions n ± (m) will be computed below. The critical behavior depends 
on the stability of the fixed points. At fixed m, for the physically relevant case m > 1, the 
e-expansion predicts four regimes: 

1) For n > n + (m), there are four fixed points, and the chiral one is the stable one. 

2) For n~(m) < n < n + (m), only the Gaussian and the Heisenberg 0(n x m)-symmetric 
fixed points are present, and none of them is stable. 

3) For n#(m) < n < n~(m), there are again four fixed points, and the chiral one is the 
stable one. For small e, the chiral fixed point has v < for m < 7 and v > for m > 7. 

4) For n < ra^(m), there are again four fixed points, and the Heisenberg 0(m x n)- 
symmetric one is the stable one. 

The antichiral fixed point is Gaussian for m — > 1 and m = —2 (or, equivalently for n — > 1 
and n = —2). Indeed, for m — > 1, u* — > for the antichiral fixed point, so that, from 
Eqs. (|3.10| ) and ( |3.11|) , we obtain r\ = and v = 1/2 at order e 3 . For m = — 2 we obtain 



u* = 3v*/2 and again rj = and v = 1/2 at order e 3 . We conjecture that this holds to all 
orders in e, and in the next Section we will provide a large-n interpretation of these results. 
The general behavior for n and m is better understood from Fig. [I]. In particular, the two 
functions n ± (m) are nothing but the two different branches of the curve that separates the 
region in which no fixed point is stable from the region in which the chiral fixed point is 
the stable one. Note that the boundaries of the different regions are symmetric under the 
exchange (n,m). Because of this symmetry it is more natural to consider the behavior in 
the variables 

S = m + n, A = m — n. (3-12) 

At fixed A there are then three regions: 



6 



20 



10 



10 



20 



Figure 1: The fixed-point structure in the (m,n) plane for d — 4. The solid line repre- 
sents the curves n ± (m). The dashed line shows n#(m) and the symmetric curve obtained 
interchanging n and m. 

1) For E > E + only the Gaussian and the Heisenberg 0(m x n)-symmetric fixed points 
are present and none is stable; 

2) For E# < E < E + there are four fixed point and the chiral one is stable; 

3) For E < E.H-, there are four fixed points and the 0(m x n)-symmetric one is the stable 
one. 

Using the above results, we can compute the e expansion of n ± (m) and n H {m). For 



n ± (m) we expand 
Then, by requiring 

and 

we obtain 



n^irn) = + nfe + n^e 2 + 0(e 3 ). 



det 



u* ,v*] n^~) = 0, 



± 



n, 



d(u, v) 
5m + 2 ± 2s, 

-5m - 2 =p —(25m 2 + 22m - 32), 
^i ± JL4 + i^C(3)±^C(3). 



24 Qx 32 sQx 8Q 2 



48 Q 2 



Here 



s = ^6(m - l)(m + 2), 



(3.13) 
(3.14) 
(3.15) 

(3.16) 
(3.17) 

(3.18) 



Ri 
R2 

Qi 
Q2 



-33024 + 18880m + 45444m 2 + 9288m 3 



1883m 4 - 417m 5 + 21m 6 + 4m 7 , 



-253952 - 160256m + 176192m 2 + 139240m 3 + 7756m 4 - 5854m 5 

-389m 6 + 58m 7 + 5m 8 , 
1632 + 1184m - 1376m 2 - 426m 3 + 31m 4 + 8m 5 , 
6176 - 2960m - 1230m 2 + 73m 3 + 20m 4 , 
(m + 8) 2 (m + 2)(m - l)(m - 7) 2 , 
(m + 8)(m + 2)(m - l)(m - 7). 



(3.19) 



For m = 2 this expression is in agreement with that given in Ref. |18fl . The expression for 
nf are singular for m = 7. However, this is not the case for nf, and indeed, by taking the 
limit we obtain 



23871617 5487 



For n# (m) we have 



1 r 



9331200 



320 



C(3). 



n H lm 



m 



4-2e + ^(6C(3)-l)e 2 + 0(e 3 



(3.20) 



(3.21) 



which is a trivial generalization of the result of Ref. [18|]. The calculation of the functions 
E + (A) and E#(A) follows the same lines. In particular, 



E+(A) = -l + h + ± (5A 2 - 24 - 2s) + 0(e 2 



(3.22) 



where s = ^6(A 2 + 18). 

From our calculation of the RG functions we can derive the fixed points of the theory. 
We expand 

M ie + M 2 e 2 + M 3 e 3 + 0(e 4 ), 



u 

V 



v 1 e + v 2 e 2 + v 3 e 3 + 0(e 4 ). 



Following Ref. [20] we define 



B 



-1 



Rmm. 



(mn + 8)(m + n - 8) 2 + 24(m - l)(n - l)(m + n - 2), 
mn(m + n) — 10(m + n) + Amn — 4, 
(m + n - 8) 2 - 12(m - l)(n - 1). 



(3.23) 
(3.24) 



We then easily find: 



--I( m + n-8)S n 



6R l/2 



(3.25) 



where we indicate by (+) the stable chiral fixed point and by (— ) the unstable antichiral one. 
In order to compute uf and vf, it is convenient to define two additional auxiliary functions: 

(377m + 14) , , /ll ^ 13 5 



Si 

So 



12 



] -ul + (m - l)(n - l)vi Q^uj - T^um + ^vj 



(3.26) 



(5mn + 82) 9 5mra — ll(m + n) + 53 
-w 2 ?Ji + 



36 



18 



-UiV l - 



13mn — 35(m + n) + 99 
72 
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Then, the 0(e 2 ) coefficients of the fixed points are given by 



uf = ±6 



1 _ 2u f- ^±f^vf) Sf - ("^K"- 1 ) (uf - v fj s. 



Ul lLmn 

2vfSf + (l - ^uf + ("-W"- 1 ^ *) S± 

U-y Ihn 



uf = ±6 i —-75 ^ . (3.27) 



The expressions for uf, vf are particularly cumbersome and they will not be reported here. 

Once the fixed points are determined, the critical exponents are computed by expanding 
in power of e the exponent series computed at the fixed point. Such a computation gives us 
the exponents only for 



> 5m + 2 + 2^6(7/1 + 2) (m-1), (3.28) 



n 
or 

n < 5m + 2 - 2^6(m + 2)(m - 1). (3.29) 

Indeed, if these bounds are not satisfied the fixed points are complex and therefore also the 
exponent series. In order to obtain series for the exponents in all the relevant domain we 
can perform the following trick. For n > n + , which is the case of physical interest, we set 
n = n + (m,e) + An and reexpand all series in powers of e keeping An fixed. In particular, 
for An = we obtain the critical exponents for n = n + . In such a case, for m = 2 we obtain: 

r] = -^e 2 + T^e 3 , 

1/ 2 V 50 50 / V 15000 15000 1000 SV ; 250 SV 7 V ; 

4 The Landau-Ginzburg- Wilson theory in the large-n 
limit 



In this Section, we study the large n-behavior of the LGW theory ( |2.2j ) at fixed m. The 
starting point is the general Hamiltonian ( |2.8|) with A?f = 0. In the high-temperature phase 
the symmetry is unbroken and thus the relevant saddle point is given by 

(S) = a, (T aP ) = 0. (4.1) 

Correspondingly, we obtain the gap equation: 

d d p 1 6<7 . 

(4.2) 



(4.3) 



(27r) d p 2 + M 2 nm 

where M 2 = r + w a. For the a-model ( |2.3| ) we obtain analogously^ 

d d p 1 771 
(2tt)<V + M 2 ~ H' 



6 This result can be obtained by using Eq. ( |1.2| ) and by taking the limit considered before Eq. (pT3|). 
Notice that, in order to keep M 2 finite in the limit, a must converge to mrji as vq — > 00. 
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(1) (2) 



(3) (4) 



Figure 2: The graphs appearing in second Legendre transform. The continuous line rep- 
resents the dressed propagator of the field, while the dashed line indicates the dressed 
propagator of an auxiliary field. 



From the gap equation we can obtain the scaling of the mass and thus the exponent v. For 
wq 7^ 0, proceeding as in the case of the ordinary 0(n) model, we obtain for 2 < d < 4 

v = — *— . (4.4) 
d-2 K J 

However, for u> = 0, we obtain simply M 2 = r , indicating 

v = I (4-5) 

for all values of the dimension d. 

Within the large-n limit we can recover the critical behavior of the theory at all fixed 
points. For generic vq and Uo, satisfying w > 0, vq > we obtain the critical chiral behavior. 
The standard Heisenberg behavior is obtained by setting Vq — 0, while the antichiral critical 
behavior is obtained at the stability boundary, i.e. by setting w = 0. It is easy to see the 
different types of excitations that appear in these cases: at the chiral fixed point both the 
scalar and the tensor degrees of freedom propagate, while at the Heisenberg and antichiral 
fixed point one observes only the scalar and the tensor degrees of freedom respectively. Note 
that, as a consequence of Eqs. (|4.4j) and ( |4.5|) , the Heisenberg and the chiral point have the 
same exponents for n = oo, and that they differ from those of the antichiral point which 
shows mean-field behavior in all dimensions. 

In order to perform the calculation we heavily relied upon the results obtained by Vasil'ev 
et al. [p2"|,p3|] , who studied the models corresponding to m = 1 with a method which lends 
itself to a reasonably simple extension, appropriate to the case we are investigating. In order 
to make our presentation self-contained, we must briefly review the essentials of the method. 

One first considers the second Legendre transform with respect to the field and the two- 



point function f24| . We indicate by _D" a,/3b (p), Ds(p), and D^ nS {p) the dressed propagators 
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of the field <p aa and of the auxiliary fields 5 and T a/S . Here a, /3, 7, and 5 go from 1 torn, 
while a and b go from 1 to n. It is useful to factorize the group dependence and to introduce 
scalar propagators 

D;^\p) = S^S ab D^p), (4.6) 

Df' yS (p) = ±(V^ 5 + 5 aS 5^ - —5 aP 5^ s ) D T (p). (4.7) 
2 V m J 

Also we can reabsorb the coupling vq and wq in the fields. Using the same notations of 
Ref. [22], we have for the second Legendre transform 

r = fir log^ + iTr lo gj D 5 + |Tr lo gj D T + f( 7l (5)+7i(T)) + 

\nm^ 2 {S) + \n(m - l)(m + 2)7 2 (T) + -nrwy^SS) + |n(m - l)(m + 2)7 3 (T5) + 
" 'm - l)(m 2 - 4) 73 (TT) + ±n 2 m 2 lA {SSS) + 



32m 

±n 2 (m - l)(m + 2)( 74 (55T) + 74 (5TT)) + — (m - l)(m 2 - 4)(m + 4) 74 (TTT) + . . . 

(4.8) 

In this equation, 71, . . ., 74 are the graphs reported in Fig. and the letters in parentheses 
indicate which auxiliary fields are propagating in the graph. In these graphs one should use 
the scalar dressed propagators and Dt while each vertex is trivially one. In the equation 
we have of course reported only those graphs that are relevant for the computation of the 
critical indices at order 1/n 2 . The group-theoretical factors in Eq. (|4.ij|) have been obtained 



by using Eqs. ( |4.6|) and ( |4.7| ) and keeping into account that the vertex (jxj/T has the form: 

(paa^pbT"* 6 -> \5 ab U^5 ps + rt^ 7 - ^ Q/3 <T 5 ) . (4.9) 

Eq. (|4.8|) is completely general and can be used in the computation of the critical indices 
for all fixed points: while we should keep into account all terms for the chiral fixed point, we 
should set T = and 5 = for the Heisenberg and the antichiral fixed points respectively. 

From Eq. (|4.8p we can derive the skeleton Dyson equations for the dressed propagators. 
It is enough to compute the variation of Y with respect to the dressed propagators. We 
obtain for the field the equation 

D^-A + u + 9l (S) + ^-(m - l)(m + 2)g x {T) + 

S 2 (55) + -(m - l)(m + 2)g 2 (ST) + -L(m - l)(m 2 - A)g 2 (TT) 
77 

+mng 3 (SS) + — (m - l)(m + 2)(^ 3 (5T5) + 2^55) + g 3 (TST) + 2g 3 (TTS)) 
—{m - l)(m 2 - 4)(m + A)g 3 (TTT) + . . . = 0, (4.10) 
while for the auxiliary fields we have 

TlTYl Tl 

D s + °s—9i + nmg 5 (S) + -(m + 2)(m - l)gs(T) + 

22 2 
n m n 



9e(SS) + — (m + 2)(m - 1)^ 6 (TT) + . . . = 0; (4.11) 



2 / 2 
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(4) (5) (6) 



Figure 3: The graphs appearing in the Dyson equations. The continuous line represents the 
dressed propagator of the field, while the dashed line indicates the dressed propagator of 
an auxiliary field. 

Dt 1 + c t + -g 4 + -g 5 (S) + — (m - 2)g 5 (T) + 

Wge(ST) + — (m + 2)(m - A)g 6 (TT) + . . . = 0. (4.12) 
8m 

Here c$ and Ct are two constants, u a momentum-independent contribution due to the 
tadpoles, A = p 2 + M 2 , and ^ are the graphs reported in Fig. |3|. As before, in parentheses 
we report which auxiliary fields are propagating. Each line is associated to a scalar dressed 
propagator, while the vertices are one. 



The critical exponents are determined following closely the method of Ref. |[22j|. As in 



Ref. [^2j we introduce two auxiliary functions: 



<x)= K * \ p(x)= \ \\ 4.13 
1 [x) Ti a a[x) 

i{ x j y) = t-— ( — dr^- ( 4 - 14 ) 

a[x)a[x — 2J 

Note the trivial symmetry of the function q(x, y) which will play a role below: 

q(x,y)=q(x,*-y). (4.15) 
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The calculation of the 1/n correction starts from assuming for x — > the following behavior 
of the dressed propagators: 



Dx(x) 



A 



x 



[l + B x x 2X + 



(4.16) 



where X is any of the fields. Here ax is related to the dimension of the field. For X = <fi 
we have = d/2 — 1 + rj/2. The correction term we report is the analytic one in the 
temperature and therefore v = 1/(2 A). From Eq. ( |4.16[ ) we obtain for the inverse functions: 



D x \x) 



A xX 2d~2a x 



1- B x q(a x ,X)x 2X + 



(4.17) 



Plugging these expressions in the skeleton equations and equating the corresponding terms 
we obtain six equations for the amplitudes. Such equations have nontrivial solutions only if 



as = «t = 2 - 7] = (3, 
and the following consistency equations are satisfied: 

M 

p{a^) = 2—p((3), 
n 

+ A) = 2, 

where M is a group-theoretical factor that depends on the fixed point: 

(chiral f.p.); 
(antichiral f.p.); 

(Heisenberg f.p. 



(4.18) 



(4.19) 
(4.20) 



M 



M+ = |(m + i; 



M~ 



2m 

M H = — 

m 



(m - l)(m + 2) 



(4.21) 



Note that M~ = for m — 1, —2, a result which follows from the fact that only a symmetric 
traceless tensor propagates. Eq. (|4.19| ) allows the determination of the first large-n coefficient 
appearing in the expansion of the exponent rj, 



(4.23) 

Vu = ~r(2-f)r(f-i)r(f-2)r(f + i)' (4 ' 24) 

and the dependence on the fixed point is encoded in the factor M. The quantity 7/n is the 
well-known result for the m = 1 model; among its most important properties we wish to 
mention that r)n — > both in the d — > 4 and in the d — > 2 limit. 



n 



rr 



rv' 



For r/i we obtain 
where 



771 = Mr/11, 
4r(d- 
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Equation ( f4.20|) should allow the determination of the exponent v. However, there is 
a subtle point that has been overlooked in the previous analyses. It is convenient for our 
discussion to introduce the auxiliary function 



q(2-r),\), 



(4.25) 



which corresponds to the left-hand side of Eq. (|4.20|) . Because of Eq. (|4.15| ), the function 
r(i], A) has the symmetry 



r(rj,X) = r(rj, § - A). 



(4.26) 



As a consequence, once r\ is fixed by solving Eq. (^4.19|) , we still have the possibility of finding 
two different solutions for A = l/2u. It is convenient to parametrize the two solutions by 



, d 
A ^4 ± 



-1+p 



(4.27) 



Using Eq. ( 4.20 ) and the fact that rj is of order 1/n, one can show that also p is of order 
1/n and therefore has an expansion of the form 



P 



Pi , P2 



n w 



+ . . . 



(4.28) 



The coefficient p\ is computed from Eq. ( f4.20|) using the fact that for n large 

A-d 



i (I - 1 + H I - 1 + p) = i (I - 1 + h 1 - p 

We then obtain 



d 



(i-&) +0( --*: 



Pl = Mp 



li 5 



where 



(d-l)(d-2) 

Pn = A _ d mi- 



(4.29) 



(4.30) 



(4.31) 



As we observed the consistency equations are satisfied for two independent choices of A. In 
order to associate to correct one to each fixed point we use the large-n estimates of the 
exponent v given in Eqs. ( [4.4|) and ( |4.5|) . Then we have 



d --l + -M+p n 
2 n 

1 - -M'pn 

n 

d „ 1 „ r u 
- - 1 + -M H p n 

2 n 



chiral f.p.; 
antichiral f.p.; 
Heisenberg f.p. 



(4.32) 



The corrections of order 1/n can be obtained by generalizing the arguments of Refs. p2 



23 1 - One must only pay attention to insert proper group-theoretical factors in front of the 



corresponding m = 1 contributions: they can be obtained from Eqs. ( (4.10| ), ( ^4.11| ), and 
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We introduce the following definitions: 

d 2 + 2d - 4 



where 



Also: 



(a) 

mi 



(b) 
V21 = 



(a) 
P21 



0) 
P21 



where 



V21 = (miY 



P21 = (Vn) 



2d(d - 2) 



, d- 1 



(d-2)(4 + 2d- d 2 )^ + 



32 + 8d - 30d 2 + 7d 3 
d(d - 4) 



(mi, 



(mi) 



M> = ^(d - 2) + ^ (2 - ^ J - V ( 5 - 2 J - ^(2). 



d _ d(6 — d) 



A-d 2(4 -d) 2 , 



(4.33) 
(4.34) 

(4.35) 
(4.36) 



2 / d(d-3) d(d-2) 
I 4-d 4-d 



d(d-2). , 2 



2(4 -d) 



;(Vu) 



d(d-2) f 

2(A-dy [r]ll) 

2d(d-3) 2 
+ l4^dT 



2(d - + -d(d - 3)R 1 + 2d - 8 + 
i^|* + 5d(3d- 8 ) Bl 



+ 



12 



4-d d-2 (d-2) 



(6i?i - R 2 - R\) + d 2 + d + 20 - 



12(4 -d) ^ 16 



4-d 



(4.37) 



R 2 = 1/(d - 3) - ^ (2 - fj - ^ - l) + ^'(1), 
^ = ^-3)+^^-^] -V^-l) -^(1). 



(4.38) 



Recalling the above considerations about the choice of A in conjunction with the choice of the 
critical point, we can now write down our final results for the 1/n expansion of the critical 
exponents in the 0(n) x 0(m) models, both for the chiral (stable) critical point and for the 
antichiral (unstable) one: 



r, 
rf 



m + 1 1 

2n ^ U n 2 



(m + 1) 2 HF m + 3 (a ) m 2 + 3m + 4 (6) 
- A V21 + — i — W + o ^21 



+ 



, (4.39) 



(m-l)(m + 2) 1 

-mi + — 

n z 



m-l) 2 (m + 2) 2 HF (m-l)(m 2 -4) (a) 



2mn 



(m - 1) (m 2 - 4) (m + 4) (b ) 
"I ~ o ^21 



8m 2 
2 m+1 



4m 2 
2 1 



-%i + 



4m 2 



(4.40) 



d-2 (d-2) 2 2n 



-P11 - 



(d - 2) 2 n 2 



m + 1) 2 HF m + 3 (a) 



-P21 + 



P21 
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m + 3m + 4 m 
+ 3 P21 



2 (m + 1) 2 



d-2 



4 



P 2 n 



1 1 (m - l)(m + 2) 1 1 

2 + 2 2n^n" Pn + 2^ 



m 



+ 
l) 2 



4m 2 



(m + 2) 2 

P21 



m 



l)(m 2 -4) (o) (m- l)(m 2 -4)(m + 4) (6) 
P21 + 5—5 P21 



+ 



m 



Am? 
l) 2 (m + 2) 2 



4m 2 







1 



(4.41) 



(4.42) 



The expressions for the stable fixed point at order 1/n coincide with those of Ref. |T^]. Note 
that i]~ = 0, v~ = 1/2 for m = —2, 1, in agreement with our e-expansion results. 

It is possible to expand the above large-n results in powers of e = 4 — d. The resulting 
expressions can be compared with the the e-expansion results for the LGW hamiltonian 
presented in the previous Section. We find full agreement both for the stable and the 
unstable fixed point for all m, thus confirming our identification of the large-n fixed points 
with the perturbative ones. 

For d — 3 the large-n expansions simplify to: 



Tj 



V 



4(m + 1) 16(m 2 - 7m - 26) 



3nn 2 



27n 2 7r 4 



4(m - l)(m + 2) 16(m - l)(m + 2)(m 2 



5WI 



3mn7r 2 
16(m + l) 
3mr 2 



27m 2 n*7r 



2-^4 



+ 



n 



1 /4(m 2 + 3m + 4) 64(5m 2 + 19m + 32/ 



7T^ 



27tt 4 



+ o(- 

m 



1 4(m + 2)(m-l) 

2 3mnir 2 



(m - l)(m + 2) 
27m 2 n 2 7r 4 



16(13m 2 + 4m + 28) + 27(m + 4)(m - 2)tt 2 



O 



(4.43) 
(4.44) 
(4.45) 

(4.46) 



5 The 1/n-expansion in the presence of a vector field 

It is quite instructive to extend the discussion of the previous paragraph to the more general 
case in which gauge-invariant vector degrees of freedom are allowed. This corresponds to 
studying the general Lagrangian (|2.1|) with to = 0. 

In the large-n limit one starts from Eqs. ( |2.8| ) and ( |2.9| ) with to = 0. As discussed by 
Hikami ||19|| , the gauge kinetic term is irrelevant for d < 4, as well as the T 2 and 5* 2 terms 
in Eq. ( |2.8| ). Thus, the large-n limit can be studied by keeping only into account H e g. The 
discussion of the fixed points is identical to that presented in the previous Section. If g 7^ 
a new set of fixed points appear: for generic v > 0, wo > we have the chiral-gauge fixed 
point in which all excitations (S, T a/3 , A°f) propagate; for wq = we have the antichiral- 
gauge fixed point (T a/3 , A°f), for v — the Heisenberg-gauge fixed point (S, A°f), and for 
u o = Vq = Wo = the pure gauge fixed point (A°f). 

Here, we want to compute the critical behavior for n — > 00, keeping only the lead- 
ing correction. Since the model is gauge-invariant, the large-n propagator of the field 
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A°f is not uniquely denned. Indeed, by integration over the field <f> we obtain a coupling 
\A°/f A/f M^ ,1& , where in momentum space 



(5.1) 



which is not invertible. A propagator for the field A°f is obtained by adding a gauge- 
fixing term, that introduces a longitudinal term, makes the matrix invertible, but does not 
contribute to physical quantities [ P5[ . 

The calculation is completely analogous to that performed in the previous Section. For 
the second Legendre transform, we obtain to order 0(1/ n) 



r = T{A = 0) + iTr logD A + § 7l (A) + 



nm(m 



-7 2 (A) + 



(5.2) 



where 71 and 72 correspond to the graphs reported in Fig. 0, all group-theoretical factors 
have been explicitly singled out, and T(A = 0) is the expression reported in Eq. 
Generalizing the results of Refs. fl^,|25]| we obtain then 



Pi = P11 



M + 



d 2 -l 



(m 



1) 



(5.3) 



where M is a group-theoretical factor defined in Eq. (|4.21|) (for the pure-gauge fixed point 
M = 0) and p\ is the 1/n contribution to the exponent p defined in Eq. (|4.27|) . Note that 
the result ( |5T3"1) does not depend on the gauge fixing used to define the propagator of the 
field A. In principle, one could also compute a gauge-fixing dependent exponent 77 but its 
significance is not so clear, since, because of the gauge invariance, the field (f) does not have 
a well-defined anomalous dimension. 

As a check we can compare our results with those obtained in perturbation theory for 
the gauge Hamiltonian (|2.1|) with to = 0. Hikami |19| determined the following one-loop f3 
functions in the MS scheme: 



(3 a 



-eu + 

-ev + 
-ea + 



mn + 8 9 
m + n 



m 



l)(n 



-v — uv 



—(m — l)ua + -(m — l)a 2 



6 

n 11 

12 ~ y' 



-v + 2uv 



-(m — l)va + -^( m — 2)a 2 



in 



a 



(5.4) 



where a = N^g 2 . These expressions generalize the results presented in Sec. p|. Choosing 
the a* = solution of the fixed-point equations we obtain the four critical points already 
discussed. However, if we choose the solution 



a 



11/ 

(m-2 



n 
12 



+ 0(e 2 



(5.5) 



we find another set of four critical points, corresponding to the distinct roots of a quartic 
algebraic equation. This equation cannot be solved in closed form, but it is easy to find its 
roots in the form of a series in the powers of 1/n. The relevant terms in the expansion of 
the roots are: 
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chiral- gauge fixed point: 



it* 1 2 - 10m , 

tt = - + — + (e , 

be n rr 

1 32 - 10m _ , 
- = - + - 2 + e . 

be n rr 



(5-6) 



antichiral-gauge fixed point: 

u* m — 1 
be" ~ 

■i>* 

be" 



+ 



(m- l)(16 + 88m- 10m 2 



ran 

1 12 + 32m - 10m : 



m 2 n 2 



+ 



0( 



n 



(5.7) 



Heisenberg-gauge fixed point: 

1 

m 

27(m - 2 



u 
be" 

be" 



1 27m 3 - 117m 2 + 90m -8 ^. , 

+ a"2 + 6 ' 

mn m z n z 



pure-gauge fixed point: 



be" 

bT 



+ 0(e 



27(m 



(5.8) 



27(m - 2) 



+ 0(e) 



0(e). 



(5.9) 



Thus, the gauge model has in general 8 fixed points and, at least for large n, the chiral-gauge 
fixed point is the stable one.[] Substituting these expressions into the relationship [19| 



mn + 2 . 
2 — — u 



(m — l)(n — 1) # 3(m — 1) 



-a* + 0(e 2 



(5.10) 



b b 4 

we find al/n expanded form of the O(e) contribution to the critical exponent v for each of 
the four solutions: 



1 

v 



2 - e + (48m - 42)- + 0(e 2 , n~ 2 ) 
n 

2 _b(m-l)(8m + l)e 



m 



chiral-gauge f.p.; 
antichiral-gauge f.p.; 



2-e + 



45m 2 — 45m + b e 



(5.11) 



rn 



n 

2 ^-2\ 



+ 0(e 2 ,n 2 ) Heisenberg-gauge f.p.; 



2-45(m- l)- + 0(e 2 ,n 



n 



pure-gauge f.p. 



It is then a matter of trivial algebra to verify that these expressions are in full agreement with 
the e-expansion of the four solutions discussed above in the context of the 1/n expansion, 
which explains the names we have given to each fixed point. Again, we think it is important 
to notice that the conformal bootstrap approach can naturally accomodate for the expansion 
of all solutions, not only the stable ones. 

7 This is not true for generic n and m. In order to obtain the general fixed-point structure, one should 
generalize the analysis performed in Sec. |^. 
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6 The e-expansion of 0(n) x 0(m) nonlinear a models 

The LGW Hamiltonian is the natural tool for the study of the critical behavior of systems 
near the upper critical dimension d = 4. If one is interested in the critical behavior near 
the lower critical dimension, one can still use perturbation theory, applied however to the 
nonlinear a models (NLer). The degrees of freedom of the NLcr models should correspond to 
the interacting Goldstone modes of the system, while the effect of the massive modes is only 
taken into account in the form of constraints for the massless fields. In this context, it is 
possible to perform an expansion in powers of e = d — 2. In the present paper we extend the 
results of Refs. [p| -|T6| to a general 0(n) x 0(m) symmetry group and to 0(e 2 ). Comparing 
with our previous 1/n expansion results we will be able to identify the nature of the fixed 
points of the NLcr. In particular, we will show explicitly that the stable fixed point of the 
generic model can be identified with the stable fixed point of the LGW theory. 



We consider the Hamiltonian ( |2.6|) . This Hamiltonian is geometric in nature, and its 
variables are best understood as generalized coordinates spanning a manifold. The cases we 
shall be interested in correspond to manifolds that are coset spaces. More specifically, we 
must study the coset space (remember that n>m) 

0{n) x 0{m) 



0(n — m) x 0(m) 

which is topologically equivalent to 



(6.i; 



°<"> (6.2) 



0{n — m) 

Associating fields it 1 with the Goldstone modes that correspond to the broken generators 

{Lie (O(n)) -Lie (0(n- to))}, (6.3) 
the Hamiltonian may be formulated in purely geometric terms, i.e. 

H = ^(tOVttW. (6.4) 

The couplings Tj = 1/rji are related to the independent entries of the tangent-space metric 

A number of important RG properties of the NLcr models have been derived in the general 
case by Friedan |2(J and specialized to the models of interest in Refs. fl5],[il|. If Rijkl and 



Ru are respectively the Riemann and Ricci tensor for the metric gu, the RG j3 functions of 
the model can be written to two-loop order as 

(3jj = s^- = -~e VlJ + R u + ±R IPQR Rj PQR + ... (6.5) 

The number of algebraically independent /3-functions coincides with the number of inde- 
pendent couplings Tj. Therefore we should consider two /^-functions associated with T\ and 
T2. The fixed points are determined from the equations 

A(T*,T*)=0, (6.6) 
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that can be perturbatively solved in powers of e. 

The evaluation of the two-loop (3 functions for arbitrary m and n requires no special 
skills, but it takes some time and effort in view of the many computational steps involved. 
Without belaboring on the intermediate steps, we report here our final results: 



Pi = S 
^2 = S 



dTi 

ds 

m 

ds 



-e7i + 
-eT 2 + 



n-2 X 



T 2 + A + BX + CX 2 T? + 0(T?) 



(6.7) 



m — 2 n — m <- 
+ — X' 



T 2 + 



DX 2 + EX 3 + FX + G 



2 2 

where X is shorthand for the ratio Ti/T^ and we have defined the coefficients: 

3 



A(n, m 
B(n, m 
C(n, m 
D(n, m 
E(n, m 
F(n, m 
G(n, m 



2m(n — m) — n H — (m — l)(m — 2), 
— {m — 1) 



Ti + 0(T?), 



(6.1 



^(n-m) + ^(m-2) 



3 . . m 
— [n — m) H 



= — (n — m)(m — 2), 
3 

= — — (n — m){m — 2), 
= — (n — m)(2m — 3), 
= i(m-2) 2 . 



Our results were submitted to a number of basic consistency checks: 

i) when m = 1 there is no rj2 coupling, and Pi(Ti) reduces to the well-known j3 function 
for the vector a-model defined on the coset space 0(n)/0(n — 1). 

ii) when m = 2 our expressions reduce to those of Ref . jnj . 

iii) when m = n there is no rjx coupling ,0 and the model reduces to a standard 0(n) xO(n) 
principal chiral model. One may verify that ^(Ti) is directly related to the known (3 function 



of these models 27 



One may also consider the "gauge" limit r/2 = 0, which was studied by Hikami ||19| : the 
identification with Hikami's coupling is 771 — > \ jt. One must however recognize that the limit 
is singular, and as a consequence the function f3l au9e (t) is not obtainable from our expressions 
by setting X = (with the notable exception of m = 2 models). If we assume 772 = from 
the very beginning of our calculation, the result is 



Pf au9e {t) = -it + {n- 2)t 2 + [2m(n - m) - n]t 3 + 0(t 4 ), 



(6.9) 



consistent with that reported in Ref. ]I9 . 



8 li n — m the vectors e Q are an orthogonal basis in R m and therefore satisfy the completeness relation 
7"^( v ■ e ")(w • e a ) = v ■ w for all vectors v, w. Then, it is a simple matter to show that the second term in 



Eq. (2^6) is one half of the first one. 
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The /^-functions ( |6.7|) are the starting point for the perturbative evaluation of the critical 
points and exponents to 0(e 2 ). A consistent Ansatz for the simultaneous solutions of the 
equations /3j(T*,T 2 *) = is the following: 



T; = tie + t 2 e 2 + 0(e 3 ), 
X* = X + X 1 i + O(i 2 ). 

It is straightforward to obtain the following algebraic equations for t\ and X 



1 m 
-- + (n-2)-- 

1 m-2 

" A » + — + 



2 

n — m 



Xn = 0, 



-x' = 0. 



(6.10) 
(6.11) 



(6.12) 
(6.13) 



They are trivially solved by 



X" 



4 



n 



2 ± yJ(n-2) 2 - (n-l)(m-2) 



n — 1 



m — l 



n 



-X*. 



Iterating the procedure we may also obtain 



± (D - A)(Xq ) 2 + (E- B)((X±Y + (F — C){X±Y + G 



±^3 



~±\4 



X±[(n - 2) - (n - l)Xt] (n 



m 



A + BXq + C(X d 



±\2 



n 



m—l x^± 



Til - lv t 
1 



L Y- 



n 



m—l 



2 ' 



(6.14) 

(6.15) 
(6.16) 



This analysis shows the existence of a couple of nontrivial fixed points of the RG equations. 
However, it is evident from Eq. ( 6.14 ) that such a pair of solutions does not exist for all m 



and n: for some values X$ is indeed complex. Repeating the analysis we performed in Sec. 
|3] for the e-expansion, we see that these two fixed points exist only for n > n + and n < n~ , 
where 



n 



m + 2 ± \Jm 2 - 4 1 
2 2' 



/m 



m 2 + 4 ± m\/ 



??? 



2(1± 



4) (m ± \/m 2 - 4) 



e + 0(e 2 ). (6.17) 



Note that for e small, we have m < h + < m + 1 and n~ < m. Thus, since n> m, all models 
with integer n > m + 1 have a a pair of nontrivial fixed points, at least for e small. Beside 
these two fixed points, there is also a fixed point for T\ — and T 2 = T 2 * that belongs to 
the universality class of the 0(m) x 0(m) principal chiral model. Such a fixed point always 
exists perturbatively for m > 2 and in particular is the only resent for n = m. 

As one may easily notice, the expansion ( p.l7| ) is singular when m — 2. This is related 
to the following peculiar feature of m = 2 models: for any value of n the (unstable) fixed 
point corresponds to the solution X* = 0, and as a consequence we observe its coalescence 
with the "gauge" fixed point obtained by setting r] 2 = 0. This phenomenon does not happen 
for m > 2. In this case, the gauge fixed point and the antichiral fixed point are distinct. 
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The exponent r/ is easily computed. To two-loop order we have 
rj = -e + (n - m)T* + (m - 1)T* + 0(e 3 ). 
Substituting the expression of the fixed point, we obtain 



Tj = (n — m){tie + t 2 e ) + (m— 1) 



X 



Ufx ' 
*o 2 . 



e 2 



-e + 0(e 



;3\ 



(6.18) 



(6.19) 



We can expand the 77 at the stable critical point in powers of 1/n, obtaining 



m + 1 3m 2 + 7m + 6 
+ ^ 



2n 



57T 



m+1 3(m + l) 2 



2n 



+ 



4n 2 



+ - 



TV 



(6.20) 

If we compare such expression with the e-expansion of r] + as obtained from the large-n 
expansion of Sec. |], Eq. ( |4.39| ), we find complete agreement, confirming the identification 
of the two fixed points. 

In NLcr models the evaluation of stability goes together with the evaluation of the critical 
exponent is, since both are related to the eigenvalues of the derivative matrix 



(T*,T*) 



(6.21) 



More precisely, stability requires that the above matrix possesses only one positive eigenvalue 
A + = is^ 1 . The presence of two positive eigenvalues signals the instability of the fixed point. 
It is possible to evaluate the above-mentioned eigenvalues in the context of the e expansion, 
obtaining: 



A_ 
A. 



where 



e - u 2 {n,m)e 2 + 0(e 3 ), 
Ai(n, m)e + X 2 (n, m)i 2 + 0(e 3 ), 

d\\ 2 a 22 \ + 0111^222 — O122O2H — 01210212 

is 2 {n,m) = ■ . 

Sill + «221 

\ x {n,m) = 1 - a m - a 22 i, 

A 2 (n, m) = is 2 (n, m) - a U2 - 0222- 



(6.22) 
(6.23) 



Here we defined 



m — 1 



am 
CZ211 

a 112 = 
^122 = 



~ tlX() — —^121^0) 

-{n - mjtxXl = -a 22 iX t 
m — 1 



(6.24) 



(6.25) 



2 

m — 1 



+ t 2 X ) - tftA + 2BX + 3CX { 



0)1 



t 2 + tl(B + 2CX ) 



(n - m)(2tiX X! + t 2 X 2 ) - i|(2DI„ + 3£X 2 + 4FX 3 ), 



«212 = 

a 222 = in- mjihXt + t 2 X ) +tj(D + 2EX + 3FX 2 - G/X 2 ). 
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The 1/n expansion of v evaluated at the stable fixed point coincides with the e expansion of 
u + obtained in Sec. |J The result of the expansion is: 



Notice that the coalescence value h can be easily determined within the e expansion by 
imposing the condition 

\_[n t ,m]=0. (6.27) 

While the stable fixed point is identified with the chiral fixed point of the 0(n) x 0(m) LGW 
model, the unstable one is unrelated to those of the LGW model. In order to understand its 
nature, it is again useful to consider the large-n limit. From Eq. ( |6.14 ) one observes that 
Xq ~l/nasn^ oo. Therefore, the fixed point survives in the large-n limit if we scale the 
coupling constants as T\ = 0(1/ n), T 2 = 0(1). Then, for large n the (3 functions decouple. 
Moreover, while f3\(Ti) gets no contributions beyond one loop (as usual in vector models), 
$2(^2) turns into the j3 function of an 0(m) x 0(m) principal chiral model [27]. Therefore, 
in this case the pattern of spontaneous symmetry breaking is highly nontrivial, even in the 
strict n — > 00 limit. This can be understood from Eqs. ( |2.£| ) and ( |2.iJ| ). In the large-n limit 
the relevant Hamiltonian is 

H = H cS + ^A 2 . (6.28) 

Now, the field (ft couples only to the gauge-invariant degrees of freedom, and thus at the 
saddle point the field Aff is a pure gauge transformation, i.e. 

Af = (0-%0) a \ (6.29) 

where O is an 0(m) matrix. Thus, for n —>■ 00, the Hamiltonian can be rewritten as the 
sum of two terms: 

H = \ £ K ■ (<9A + M 2 )<t> a + |Tr d,0~%0, (6.30) 

where, as in Sec. f|, M 2 = r + w a. Thus, the unstable fixed point is directly related to the 
nontrivial fixed point of the principal 0(m) x 0(m) chiral model. 

Finally, we consider the gauge limit 771 = 1/t, i] 2 = 0. We find a nontrivial fixed point: 

e + -, — — J -i 2 + 0(i 3 ). (6.31) 



n - 2 (n - 2) 

Correspondingly we obtain: 



1 r,,, *n ~ 2m(n — m) — n r, ^ . n. 
v- 1 = p(t*) = e + \ -4 e 2 + 0(e 3 ). (6.32) 



(n - 2) 2 

The 1/n expansion of the above result gives 



v- 1 = e + ^-^e 2 + 0(f), (6.33) 



23 



and one may easily check that it agrees with the e expansion of the results found in Sec. |5] 
for the stable fixed point of the gauge model, see Eq. fl5.3|), with M = M + . 



The behavior of NLa models in the case m = n — 1 is worth a special discussion [^8 
In this case one may naturally define a new n-component field e n such that e„ • e n = 1 and 
e Q ■ e n = 0, and one can show that 



i 2=K n i 



When 7]i — r)2, this is the Hamiltonian of an 0(n) x 0(n) principal chiral model. The stable 
fixed point is characterized by the property that X* = 1 and one finds: 

h = h = l —. (6.35) 

n — 2 n — 2 

Direct substitution shows that 

z/- 1 = e + ^e 2 + 0(e 3 ), (6.36) 

As one may easily check, these exponents coincide with those obtained in the case m = n. 
Therefore, as shown in Ref. [p8|l , the symmetry of the 0(n) xO(n-l) model is dynamically 
promoted to 0(n) x 0(ri) at the stable fixed point. This property is certainly true for 
sufficiently small d > 2, but at this level of analysis it is impossible to establish the maximum 
dimension d c for which a stable critical point possessing the enlarged symmetry can be found. 



7 Conclusions 

At this stage of our analysis, we can draw quite general conclusions on the general fixed-point 
structure of the models with 0(n) x 0(m) for all dimensions 2 < d < 4. By comparing the 
e-expansion results near two dimensions, the e-expansion results near four dimensions and 
the large-n results, we have been able to identify the nature of all (stable and unstable) fixed 
points of these models. In particular, the LGW stable fixed point coincides with the stable 
one of the NLcr model. We thus quantitatively confirm one of the conclusions of Ref. \ TE , IE \: 



above n(m, d) in the (m, d) plane a second-order phase transition occurs, which is described 
by the common fixed point of the NLa and LGW models. 

The unstable fixed points, which give rise to different types of tricritical behavior and 
crossover phenomena, are instead unrelated and correspond to systems with completely 
different types of excitations. 

The correspondence we have found holds only for sufficiently large values of n, i.e. for 
n > n(m,d), which is the region of analyticity of the large-n expansion. Since the 1/n 
expansion commutes with the e = 4 — d expansion of the LGW Hamiltonian, one may expect 
n(m,d) to coincide with n + (m,d) in a neighborhood of d = 4; n + (m,d) might in turn be 
evaluated within the 1/n expansion by solving the coalescence equation: 

u + {m,n + ,d) = v~ (m,n + ,d). (7.1) 
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m 



Figure 4: Sketch of the large-n analyticity boundary as a function of the dimension d. 



The estimate obtained from the lowest order approximation for z/ ± is:[] 

n + (m,d)^2(m+l--)^l. (7.2) 
V mj A- d K ' 

This expression shows the correct qualitative behavior for all 2 < d < 4 and a rough 
quantitative agreement. It is possible to improve the approximation by including the 1/n 2 
correction in Eq. ( |7.1| ). For m = 2 and d — 3, it predicts n + (2,3) ~ 5.3, in substantial 
agreement with the results obtained by using the ERG approach M, [| , the perturbative 
expansion in fixed dimension flIU| , and, as we shall show below, the e-expansion. 

We want now to understand the behavior of n(m, d) near two dimensions. Near two 
dimensions, using the NLa model results we know that the (LGW and NLa) stable fixed 
point exists only for n > n + (m, d), so that in this case n(m, d) = n + (m, d). Thus, for generic 
values of d we conjecture 

_/ ,\ | n + (m, d) for dJm) < d < 4 oN 
n(m,d) = {^ + ) ' o\a7aT\ 7 -3 

y n + (m, d) for 2 < d < d c {m), 

where d c (m) is a critical dimension that we cannot determine with our means. Of course, 
this expression is valid for n > m. The symmetry under exchange of n and m, implies the 
existence of a similar boundary curve in the region n < m, obtained by interchanging n and 
m. A sketch of n(m, d) is reported in Fig. ||. 

Now, let us discuss the behavior of the LGW fixed points for d —>■ 2. Since the LGW 
stable fixed point is equivalent to that of the NLcx model, and, for all n > 2, m > 2 except 
n = 2, m = 2, the NLcx model is asymptotically free, we expect v + = oo, a conclusion that is 
confirmed by the large-n expression ( 4.41 ). On the other hand, for d = 2, Eq. ( ^4.42|) predicts 



v =1/2 without 1/n and 1/n 2 corrections. It is thus natural to conjecture that v~ = 1/2 
for all n > 2 and m > 2, i.e. that the LGW antichiral fixed point is a Gaussian fixed point. 
The case m = 2, n = 2 needs a special discussion. Using the fact that the 0(2) x 0(2) LGW 



9 This expression has been obtained by solving exactly the equation l/V + (m, n + , d) = \jv (m,n + ,d), 
where \/v + (m, n + , d) and l/v~(m, n + , d) are expanded to order to 1/n. 
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model is equivalent to the so-called mn model with m = n = 2 [|T^,^9j, one can show that 
in the v < region a stable fixed point exists for all values of d [1110. Finally, for d — > 2, 
using the cr-model results of the previous Section, one finds that it becomes Gaussian. 

We want now to use the knowledge of n(m, 2) in order to obtain some informations on 
n(m, 3). For this purpose we will make two hypotheses: first we will assume d c (m) < 3, so 
that n(m, 3) = n + (m, 3); second, we will assume n(m, d) to be sufficiently smooth in d at m 



fixed, so that we can use the interpolation method of Ref. [|3T]]. Such a method has provided 



very precise estimates of critical quantities (see, e.g., Refs. [pl] , |32|] ). 
Let us first consider the case m = 2. We start from | fL8|1 



n(2,4 



12 + 4^- [ 12 



+ 



137 91 r /13 47 rx w , 

+ V6 + — + — V6 C 3 

150 300 V 5 60 1 



e 2 + 0(e 3 ) 



21.80 - 23.43e + 7.09e 2 + 0(e 3 



(7.4) 



Following Ref. [31], we rewrite this equation in the following form 



n(2, 4 - e) = 2 + (2 - e) (9.90 - 6.67e + 0.16e 2 ) + 0(e 3 ). 



(7.5) 

Note that the new perturbative series is much better behaved than the original one, the 
coefficients of the series decreasing rapidly. Setting e = 1, we obtain an estimate for n(2, 3): 



fi(2,3)« 5.3(2), 



(7.6) 



where the "error" indicates how the estimate varies from two loops to three loops. It should 
not be taken seriously; it should only provide an order of magnitude for the precision of the 
results. The estimate (|7.6|) is in good agreeement with the determinations of Ref. [§, [10]] : 
n(2, 3) ~ 5 (Ref. ||), ~ 6 (Ref. flOfl ). We can try to estimate the exponents for n = n(2, 3), 
by using Eq. (|3.30| ). The coefficients decrease steadily with e and thus we can simply set 
e = l, obtaining 

n w 0.038, v w 0.63. (7.7) 

It is also interesting to compute the exponents for n = 6, m = 2, and d = 3, in order 
to make a numerical comparison with the results of Ref. |3^] who found v = 0.700(11), 
7 = 1.383(36), and the ERG results of Ref. § who found v w 0.707, 7 w 1.377. If we use 
our 0(n~ 2 ) expansions for the critical exponents, we obtain 7 pa 1.22 and v ~ 0.63. We 
can also use the e-expansion, by using the method explained at the end of Sec. B]. In this 
case, we must fix An = 6 — n(2,3). Conservatively, we have < An < 1. Then, from 
the perturbative series we estimate v ~ 0.63 - 0.64, 7 pa 1.24 - 1.26, which is rather close 
to the large-n result, and somewhat lower than the numerical results of Refs. |9|,[33|]. It is 
also worth mentioning that for n = 6 the fixed-dimension field-theoretical approach do not 



10 The mn model with m = 2 describes n XY models coupled by an 0(n)-symmetric interaction. Using 
essentially nonperturbative arguments (see e.g. Ref. |2{|) related to the specific- heat exponent of the XY 
universality class, one can argue that for d < d c with d c > 3 there is a stable fixed point belonging to the 
XY universality class, while for d > d c there is a stable fixed point with the tetragonal symmetry. This fact 
has been recently confirmed by high-order field-theoretical calculations in three dimensions |0,pQ]. 
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find fixed points that are sufficiently stable with respect to the order of the expansion up to 



six loops |L(J. We believe that these apparent discrepancies among the various approaches 
deserve further investigation. 

Our expressions may also be employed in order to compute the critical dimensionality d c 



such that, for d < d c the 0(3) x 0(2) has a nontrivial fixed point with symmetry 0(4) |15,16 



This corresponds to solving the equation n(2, d c ) = 3. If we use Eq. (|7.5|), we find d c ~ 2.71 



We may compare our result to those obtained in the ERG approach: d c = 2.83 (Ref. @), 
2.87 (Ref. ||). They are in substantial agreement with our result, when allowing for the 
systematic errors of both approaches. It should also be noticed that our interpolation ( |7.5| ) 
is also in very good agreement with the ERG results of Refs. [^,|| for all e [|34| . 



These analyses can be repeated for larger values of m. Since only n(m, 3) seems to be 
rather precisely determined, we only report the results for this quantity. For m = 3 and 4 
we have the constrained estimates 

n(3,3) « 9.1(9), (7.8) 
n(4,3) « 12(1). (7.9) 

For large m we have 

n(m,4-e) = m (9.90 - lO.lOe + 2.66e 2 + 0(e 3 , m" 1 )) 

= m + (2-e)m [4.45 - 2.83e - 0.084e 2 + 0(e 3 , m" 1 )] , (7.10) 

where, as already observed, the coefficients of the constrained series are smaller than the 
original ones. Setting e = 1, we obtain n(m, 3) ~ 2.5m. Note that the large- m approximation 
is already good at m = 4. 

It is very important to notice that, since all extrapolation techniques are adiabatic in 
their parameters, it is not possible to catch the (essentially nonperturbative) features of 
the models in the region below n. As a consequence there is no inconsistency between the 
present statements and our results ]10|] concerning 0(n) x 0(2) models for n = 2, 3 in fixed 
dimension d = 3. The fixed points we found for n = 2, 3 are certainly not analytically 
connected with the large-n and small-e criticalities discussed in this paper. 

Finally, we recall that an enlarged parameter space for the 0(n) x 0(m) symmetric 
models with critical dimension d = 4 leads to the appearance of several new, generally 
unstable, fixed points, that physically correspond to tricritical transitions and give rise to 
crossover phenomena. It is important to recognize that the conformal bootstrap approach 
to the 1/n expansion allows a consistent treatment of all these criticalities. Systems with 
0(n) x 0(m) symmetry may also possess a "gauge" criticality, which can be described by the 
appropriate 1/n expansion as well as within the e expansion of the Hamiltonian for scalar 
chromodynamics and within the e expansion of a class of gauge-invariant NLcx models. 
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